Precision calculations of the gravitational wave background spectrum from inflation 
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The spectrum of the gravitational wave background originating horn quantum fluctuations during 
inflation is calculated numerically for various inflation models over a wide range of frequencies. We 
take into account four ingredients : the scalar field dynamics during inflation making no use of 
the slow-roll approximation, the fermionic decay of the scalar field with a small coupling constant 
during the reheating process, the change of the effective number of degrees of freedom p, in the 
radiation-dominated era, and the anisotropic stress of free-streaming neutrinos. By numerically 
solving the evolution of gravitational waves during and after inflation up to the present, all of 
these effects can be examined comprehensively and accurately over a broad spectrum, even at very 
high frequencies. We find that the spectrum shows (i) a large deviation from the spectrum less 
accurate obtained by Taylor expanding around the CMB scale using the slow-roll approximation 
(ii) a characteristic frequency dependence due to the reheating effect, and (iii) damping due to the 
p* changes and the neutrino anisotropic stress. We suggest that future analysis of the gravitational 
wave background should take into consideration the fact that analytical estimates using the Taylor 
expansion overestimate the amplitude of the spectrum. 
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I. INTRODUCTION 

Inflation [l|, 0, W\ is now widely considered to be the 
most attractive scenario for describing the early evo- 
lution of the Universe. One of the reasons why in- 
flation succeeds as a standard paradigm is that it not 
only solves conundrums in the big bang theory (hori- 
zon/flatness/monopole problems) but also predicts an al- 
most scale-invariant spectrum of scalar (density Vpertur- 
bations originating from quantum fluctuations [1, H, Q . 
This theoretical prediction is consistent with large-scale 
structure and cosmic microwave background (CMB) ob- 
servations, and there have been intensive attempts to 
obtain a deeper understanding of inflation dynamics 
through analysis of these observations. Moreover, infla- 
tion also generates an almost scale-invariant spectrum 
of tensor perturbations (gravitational waves) 0, [a, Q, 
the detection of which is expected to provide a new 
source of information on the early Universe, and pro- 
vide a key to discriminate between inflation and alter- 
native models. There are ongoing efforts to detect in- 
directly this inflation-produced gravitational wave back- 
ground, such as the next-generation CMB experiments 
which measure CM B p olarization patterns induced by 
gravitational waves Il0l[l]| . Clearly, the direct detection 
[12 . [l3l . UA . [la . UM . [171 of the primordial gravitational 
wave background would provide an even more attractive 
prospect. 

The most striking property of gravitational waves is 
their weakness of interaction with matter. While this 



property makes direct detection very difficult, it may 
provide the following specific benefits. Since primordial 
gravitational waves propagate freely since the end of in- 
flation because of their weak interactions with matter, 
one may observe a remnant of inflation directly through 
observation. In addition, due to the smallness of their 
amplitude, a linear approximation is valid for gravita- 
tional waves in all scales. This makes it easy to study 
their behavior at any frequency. This is in contrast to 
density perturbations, whose evolution is more compli- 
cated, and which have now grown nonlinear at small 
scales. [Although gravitational waves are affected by the 
second order term of the primordial dens ity p erturbations 
[18J and the growing nonlinear modes [19[ through the 
anisotropic stress term, these effects appear at very low 
frequencies ^ (< 10~"'^^IIz) and have been estimated to be 
small.] Therefore, the gravitational wave background has 
the potential to allow us to extract "clean" information 
from the early Universe, and is extremely attractive as a 
tool of future observational cosmology. 

In this paper, our goal is to provide precise predictions 
for the spectrum amplitude of the gravitational wave 
background as preparation for assessing detectability in 
future experiments aiming at the direct detection of grav- 
itational waves. It is a well-known fact that the tilt and 
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^ One might expect more gravitational waves are generated on 
smaller scales due to the nonlinear structure formation since the 
structure becomes nonlinear from the smaller scales in the CDM 
model. However, it is not the case because the amplitude of the 
gravitational waves generated at high redshift decreases while 
they propagate in the expanding Universe, and becomes negligi- 
ble small. (See Ref. Il9l . for details.) 



the amplitude of the spectrum strongly depend on the in- 
flation model. Also, after inflation, changes in the Hub- 
ble expansion rate and the anisotropic stress of neutrinos 
induce characteristic features in the spectrum [20, l2]| . 
These components which affect the spectrum should be 
considered all together to estimate accurately the ampli- 
tude of the spectrum over a wide frequency range. The 
latter effects have been well investigated at frequencies 
below 10~^Hz by Watanabe and Komatsu [22, who cal- 
culate the spectrum taking into account both the change 
of the effective number of degrees of freedom and the 
neutrino anisotropic stress. In contrast, our interest lies 
particularly in the higher-frequency region of the gravi- 
tational wave background spectrum. 

The reason why we focus on the high frequency modes 
is that they are sensitive to a crucial aspect of the dy- 
namics of inflation: how inflation terminates. The high 
frequency modes are produced outside the horizon near 
the end of inflation and reenter the horizon soonest after 
inflation ends, in the period where reheating is consid- 
ered to occur. Since the amplitude of a gravitational- 
wave mode is affected by the Hubble expansion rate at 
the point when that mode crosses the horizon, the spec- 
trum of these modes depends on the particular inflation 
and reheating models chosen. This means that it may be 
possible to extract information about the epoch, which is 
still quite unclear, through direct detection of the gravi- 
tational wave background. 

In this paper, we numerically analyze the spectrum of 
the primordial gravitational wave background for several 
inflation models. It is particularly worth noting that the 
dynamics of inflation and reheating are included by nu- 
merically solving the relevant equations, which describe 
an inflaton field rolling down its potential and subse- 
quently decaying into radiation with oscillation at the 
bottom of the potential. This means that we make no use 
of the slow-roll approximation, which is no longer valid 
near the end of inflation, unlike many previous works 
[23, [2J, [25, 26]. Our numerical approach for computing 
the dynamics of inflation enables us to find the devia- 
tion from this approximation, which depends strongly on 
the inflation model. Moreover this numerical approach 
works effectively to connect the change of the Hubble 
expansion rate smoothly from the inflation phase to the 
radiation-dominated era via the reheating process, which 
cannot be accurately treated analytically. Since the Hub- 
ble expansion rate depends on the behavior of the scalar 
field and following reheating process, the amplitude of the 
high frequency modes are determined by a combination 
of both the dynamics of inflation and the reheating mech- 
anism. Detailed numerical calculations of these processes 
are therefore required in order to examine precisely the 
high frequency gravitational waves which are the main 
subject of this paper. 

In our calculation, we take into account the change of 
the effective number of degrees of freedom and the neu- 
trino anisotropic stress, both of which cause a damping 
of the spectrum. It is notable that our detailed calcu- 



lations have also revealed a numerical error due to the 
inaccurate treatment of the effect of neutrino anisotropic 
stress in previous work (Ref. Q). As a result, we obtain 
a wide range spectrum which shows interesting features 
produced in the early Universe, the inflationary phase 
and the reheating phase. We should therefore stress that 
this is the first work which calculates the amplitude of the 
gravitational wave background accurately including all 
of the factors that are definitely understood to affect it. 
This makes it possible to obtain an accurate spectral am- 
plitude and is therefore essential for accurately assessing 
the detectability of the gravitational wave background. 

The outhne of our paper is as follows. In Sec. m we be- 
gin with the evolution equations expressing gravitational 
waves in terms of metric tensor perturbations. Then we 
review the basics of the gravitational wave background 
spectrum and evaluate the basic spectrum shape analyt- 
ically. In Sec lIIIl we describe the physical phenomena 
in the early Universe that are thought to have affected 
the spectrum shape. In Sec lIVi we present the results 
of our numerical calculations of the gravitational wave 
background spectrum and provide comparisons with the- 
oretical predictions. We also show the results at different 
reheating temperatures and different inflation models. In 
Sec|Vl we conclude with a discussion of the results of 
SecHVl 



II. THE EVOLUTION OF GRAVITATIONAL 
WAVES IN THE EXPANDING UNIVERSE 

Gravitational waves in the expanding Universe can be 
described as a transverse-traceless part of the metric per- 
turbation in a Friedmann Robertson- Walker background. 
We consider the spatially flat case in which the line ele- 
ment is given by 



ds'' 



-df + a^{t){5r 



hij)dx^ dx-' , 



(1) 



where the tensor perturbations hij satisfy the transverse- 
traceless conditions, ft-00 = hoi = d'^hij = hi = 0. Apply- 
ing the metric of this space-time to the perturbed Ein- 
stein equation SGij = SnGSTij, we obtain the equation 
of motion for gravitational waves 

hij + iHhij ~ —^"^hij = 167rGHy , (2) 

where the over-dot describes the time derivative, H is 
the Hubble parameter, and Hy is the transverse-traceless 
part of the anisotropic stress of the energy-momentum 
tensor. The obvious contribution from the anisotropic 
stress term is the damping effect due to neutrino free- 
streaming, presented in Rcfs. [U, mm, [23. We take 
into account this effect numerically solving the integro- 
differential equation shown in Appendix [C) Although 
there may be other yet unknown free-streaming particles 
which cause damping at higher frequencies, we do not 
include such an uncertain component. 



We now develop the basics of the gravitational wave 
background spectrum assuming a perfect fluid 11^ — 
for simplicity. For convenience, we decompose hij into 
its Fourier modes with the comoving wave number k and 
denote the two independent polarization states as 



htj{t,x.) 



A=+,x 



(27r)3/2 



,(k)/i^(t,k)e' 



(3) 



In the second step, we have substituted the Fourier trans- 
formed form of /OgWj which is given by the 00-component 
of the stress-energy tensor as follows, 
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where the polarization tensors e^ ■' satisfies symmetric 
and transverse-traceless condition and are normalized as 
J2^ . e^Ae^jY = 26^^'. Then Eq. © becomes 






(4) 



Note that, by neglecting the anisotropic stress term, the 
contribution of other components in the Universe appears 
only through the second term. This means we can follow 
the evolution of gravitational waves with Eq. ^ simply 
by computing the Hubble expansion rate, the evolution 
of which is determined by the total energy densities of 
the background components. 

Let us briefly explain how the primordial gravitational 
waves evolve in the expanding Universe. Initially, dur- 
ing inflation, the modes of interest were deep inside the 
Hubble horizon (fc ^ aH) and fluctuating quantum- 
mechanically. During this phase, the amplitude of h^ 
is determined according to Appendix [XJ 



\h'^\' 



IGttG 

2fca2 ■ 



(5) 



As the Universe expands exponentially, the mode crosses 
outside the horizon {k < aH). At this moment, we 
assume that the quantum fluctuations become classical. 
Using the fact that the second term in Eq. (j4|) becomes 
dominant compared with the third term, one can easily 
obtain the solution 



h^ ex const. 



(6) 



After the end of inflation, the mode reenters the hori- 
zon (fc > aH) and starts to oscillate. In this phase, the 
behavior of h^ is described by the WKB solution (see 
Appendix [B| , 



/i^ oc a-le±*'=^ 



(7) 



where r is the conformal time defined by dr = dt/a{t). 

We now define a dimensionless quantity to characterize 
the strength of gravitational waves in terms of their en- 
ergy density pgw and the critical density of the Universe 
Pc = 3H^/8ttG [li], 



n 



GW 



1 dpi 



'GW 



Pc d\nk 
1 / k 
12 V^ 



fc3 



E 



h'^\'- 



(8) 



Here, we have assumed that the mode is inside the hori- 
zon and behaves like Eq. ([7]), which means that the 
conformal-time derivative of the mode can be regarded as 
the same as the spatial derivative when taking the spatial 
average, (■••). 

One may define the tensor power spectrum Vrik) in 
terms of the primordial tensor power spectrum 7^T,prim(^) 
and the tensor transfer function TT{k) as 



fc3 

7r2 



'^T{k) ^'-^Y. I'^kl' = ^T,pnm(A:)r|(fc). (10) 



Then the energy spectrum ficw can be rewritten as 



VL 



GW-Y^(^) PT,prim(fc)T|(fc). (11) 



This expression is used in previous works because it is 
useful to evaluate the spectrum by splitting the contri- 
butions to the spectrum shapes into two parts, VT,pmm{k) 
and Txik). The spectrum shape is affected by the Hub- 
ble rate both when the mode leaves the horizon during 
inflation and when the mode reenters the horizon after 
inflation. Since Vr.primik) is the spectrum at the end 
of inflation, it contains only the contribution from the 
former, whereas TT{k) represents the time evolution of 
each mode after the end of inflation, which means it con- 
tains only the contribution from latter. In most previous 
work, 7^T,prim(fc) has been evaluated analytically using 
the slow-roll approximation. Wc shall return to this in 
the next section. 

We now show two simple examples of frequency de- 
pendence of the spectrum - a radiation-dominated and a 
matter-dominated case - assuming the primordial spec- 
trum is flat, ■PT,prim(fc) oc /c". Accordiug to Eq. ©, h^ 
remains constant on superhorizon scales. This means the 
transfer function can be written in terms of the ampli- 
tude of h^ at horizon crossing, T^{k) — |/ik,oP/|^k,hcP, 
where the subscript "0" denotes the present time and 
" he" denotes the horizon-crossing time. After the modes 
reenter the horizon, h^ decreases inversely proportional 
to the scale factor, so that l/i^ol = I'^k hcK'^hc/'^o), 
then we get T^{k) — a'^c/^o- Therefore, the energy 
spectrum at the present time is rewritten as r^Gw.o = 
(fc2/12ao.ffo)^T,prim(ahc/^o)- ^y using the relation k = 
ahcH\ic, the value Ohc in f^GW.o can be converted into 
k with the proportional relation between ahc and i?hc- 
During a radiation-dominated era, the Hubble rate be- 
haves as -ff ex a~2. This leads to ahc oc k~^, then we 



obtain 



f^GW.o « k°. 



(12) 



Similarly, during a matter-dominated era, H (x a '^" 
leads to Ohc oc A;^^, then we obtain 



^GW,o oc k' 



(13) 



III. 



PHYSICS OF THE EARLY UNIVERSE AND 
ITS EFFECT ON THE SPECTRUM 



As was seen in the previous section, the evolution of 
the Hubble expansion rate is important to determine the 
shape of the gravitational wave background spectrum. 
Here, we shall describe the physics of the early Universe 
which affect the spectrum shape and provide the equa- 
tions necessary for calculating the evolution of the Hub- 
ble rate. 



A. Slow-roll inflation 

In most models of inflation, the accelerated expansion 
is driven by a scalar field </) whose equation of motion is 
given by 



3H^ +V' ^0, 



(14) 



where prime denotes derivative with respect to (j). The 
Hubble parameter is determined by the energy density of 
this scalar field, p^ — <p^ /2 + V, so that the Friedmann 
equation can be written as 



H' 



Stt 



3m% \2 



02 + y(0) 



(15) 



where mpi = l/vG is the Plank mass. For exponential 
expansion, the energy density of the field should be 
dominated by its potential energy, (j)'^/2 <C V. This con- 
dition is often characterized by the slow-roll parameters 
which are defined as 



m 



PI 



V 



XL 

Stt V 



(16) 
(17) 



As long as the so-called slow- roll conditions, e ^ 1 and 
I77I <C 1, are satisfied, inflation continues keeping the 
Hubble rate nearly constant and we can use the slow- 
roll approximation which is a useful tool to theoretically 
study the fluctuations generated during inflation. 

In the slow-roll regime, the spectra of the primordial 
scalar perturbation 7-'s,prim(^) and tensor perturbation 
7^T,priin(fc) which are generated outside the horizon can 



derived as 
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(18) 


7'T,prim(fc) = 


16 

TT 


\rnpij 


1 
k=aH 


(19) 



where H is evaluated at the time when the mode with 
wave number k exits the horizon. To parametrize the 
deviation from scale invariance of the tensor spectrum 
VT.primik), One may use the tensor spectral index ut 
and its running ax which can be expressed in terms of 
the slow-roll parameters [2^, |30| , 



nxik) 
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_ dlnVT.pnmik) 



dink 



-2e, 
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ie\ 



(20) 

(21) 



Then the power spectrum is represented in an expanded 
form in terms of ln(fc//co), 



, ^T,prim(fc) ., x, k 

^^ -n (I, \ = "^(^=0) In -r- 

PT,piim[ko) ka 



-aTiko)\n^ — 

Z fco 



(22) 

where feg is a pivot wavenumber. Also, the tensor-to- 
scalar ratio can be written in terms of e. 



VT.pnmjk) 
' S',primv"'J 



16e. 



(23) 



These expressions are used as a reference for comparison 
with our numerical results in the subsequent section. 



B. Reheating process 

During reheating process that follows inflation, the en- 
ergy of the scalar field is transferred to radiation or other 
light particles oscillating near the minimum of its poten- 
tial. Here, we consider the perturbative decay of the 
scalar field [3l|, [34I and only consider the interaction be- 
tween the scalar field and light fermions. Then the equa- 
tions for can be simply rewritten including the decay 



term as '33', '341 



4)+{3H + T)^ + V' ^0, 



(24) 



where F is the decay rate of (p particle, which depends 
on oscillation frequency at the bottom of the potential. 
Note that this expression is valid only for the case where 
the mass of the interacting field and the decay rate F is 
sufficiently smaller than the oscillation frequency of the 
(j) field. We do not deal with nonlinear processes like 
preheating [3^, [3y| , which occurs if the decay is bosonic 
and the parameters satisfy specific conditions [3J|. Such 
processes may also generate gravitational waves at high 
frequencies with a quite different spectrum shape, and 



may even exceed the amplitude of the inflation-produced 
gravitational wave background [33, S, ll^- However, 
since we aim to provide a conservative estimate of the 
spectrum, we only consider perturbative processes which 
exist regardless of the chosen reheating scenario. 

The additional set of equations to describe the reheat- 
ing process is as follows. The evolution equation for the 
radiation energy density pr can be written as 



as 



Pr + AHpr = Tp^. 



(25) 



The Friedmann equation is written in terms of the energy 
density of the field and radiation field, 



H^ 






{P<)> + Pr) 



(26) 



Therefore, using this sets of three equations, Eqs. 
([M)l . (P5l) and ([^5]) . we can simply calculate the Hubble 
rate following the evolution of both the scalar field and 
the radiation field which are interacting with each other. 
Reheating ends around T ^ H where the field be- 
gins to decay exponentially. At this stage, the Universe 
becomes radiation-dominated, so that we can define the 
reheating temperature which can be related to T [33|, 



T, 



RH 



£) ^^^' 



(27) 



where 5, (T) is the effective number of degrees of freedom. 
This value is often used to characterize the reheating pro- 
cess. 



C. Thermodynamics of the radiation era 

After the Universe becomes radiation-dominated, the 
Hubble rate basically behaves as iJ ex a~^ . In the early 
stages of the radiation-dominated era, particles in the 
Universe are relativistic and contribute to the radiation 
energy density. However, when a particle species becomes 
nonrelativistic as the Universe expands and the temper- 
ature decreases, the contribution of the particle species 
to the radiation drops exponentially and becomes negli- 
gible. This results in the change of the radiation energy 
density and the evolution of the Hubble rate, then in- 
duces a damping of the amplitude at high frequencies. 
These processes for all particle species are treated in the 
framework of thermodynamics as changes of the "effec- 
tive number of relativistic degrees of freedom," (7*(T), 
which is introduced for the convenience to evaluate the 
contribution of all relativistic species as a whole. We take 
into account this effect with a formulation based on Ref. 
[33| which is briefly described as follows. 

During the radiation era, matter particles are in ther- 
mal equilibrium and they obey Fermi-Dirac or Bose- 
Einstein statistics. In this case, the energy density and 
the entropy density at temperature T can be expressed 



p{T) ^ ^9*{T)T\ 



s{T) 



2^ 
45 



9*s{T)T^ 



(28) 
(29) 



Here g* (T) and gs * (T) can be written in the form of a 
summation of the contributions of each species as 



m-Y.9*AT)[^ 



T, 



9*s(T) ^^g^sAT) 



(30) 



(31) 



where g*^i{T) and g^s,i{T) represent the contribution of 
the ith particle species with mass rrii, degrees of freedom 
gi and temperature Ti, and are given by 



9*AT') 



15 



^{T)=9. 



15 



{u 



e"±l 



e"±l 

4 



■u^du, (32) 



du. (33) 



For this expression, we define Xi = mi/T and u — Ei/T ^ 
where Ei = -^/Ip^p -I- m|, and assume that the chemical 
potential is negligible. The sign -I- in the denominator 
denotes fermions and — denotes bosons. 

With the help of entropy conservation, sa^ — const, 

Eqs.jSH]) and ^ lead to p oc 5,57/^^0-^. This 
means the behavior of the Hubble rate in the radiation- 
dominated era is modified as, H (x gj g^,s a~^. Let 
us briefiy consider how the the changes in g^ and g^s 
affect the spectrum of the gravitational wave back- 
ground. Referring back to the discussion at the end 
of Sec. ini the energy spectrum is written as JIgw-O = 
{k'^ /12alHg)'PT,piim{a1c/al). Substituting the modified 
Hubble rate into the relation k = ahc^hc, we obtain 



1/2 -2/3, 



-1/3, 



aiic oc 5*;hc5*s,hc^ ^- This leads to flowfi « g^^^^ k°, 
where we assume g^, — g^.^, which holds before the neu- 
trino decoupling, for simplicity. Therefore, since g^ is 
large at the early stage of the Universe when the temper- 
ature is high enough for particles to be relativistic, the 
spectrum is suppressed at the higher frequencies which 
enter the horizon earlier. 

Taking into account the improved evolution of the en- 
ergy density, the Friedmann equation for the Hubble rate 
can be expressed in terms of the density parameters as 



H = Hn 



9* 
9*0 



9*3 

g*so 



nr — 



\ If] 



n^ 



(34) 



where $7^, ^m, f^nd Q\ are the density parameters for 
radiation, matter and the cosmological constant, respec- 
tively. Throughout this paper, we adopt ilrh^ = 4.31 x 



10~^ and the WMAP cosmological parameters h = 0.732, 
ilm = 0.241, and n^ = 0.759 ^]. The particles we as- 
sume are hsted in Table I of Ref. [22]. Here, we con- 
sider only particles in the standard model and do not 
include SUSY particles or any other exotic particles. In- 
cluding SUSY particles would double g* at above the TeV 

— 1/3 

scale and r2Gw,o(c< 3» he ) would decrease by a factor of 
2-1/3 -0.8 fH. 



IV. SPECTRUM OF PRIMORDIAL 
GRAVITATIONAL WAVES FROM INFLATION 

Before showing our numerical results, we briefly de- 
scribe details of the numerical calculation method. We 
first compute the background equations for the Hubble 
rate using the procedure described in Sec. Illli and then 
substitute it into the evolution equation of gravitational 
waves, Eq. (g]), presented in Sec. |TT1 Note that this equa- 
tion is obtained by neglecting the anisotropic stress term, 
which can be set to zero when neutrinos do not induce an 
anisotropic stress. We take into account the contribution 
of the neutrino anisotropic stress by using Eq. (JC1|) after 
the time of neutrino decoupling. 

For the calculation of the Hubble rate, we use Eqs. 
(IM1),(I2S1) and (gni) during inflation and reheating. Then 
we switch to Eq. (|34p at the point when the Universe 
becomes well dominated by radiation, H ex a~^. The 
e-folding number of inflation, corresponding to the pivot 
scale, is chosen in order to ensure smooth connection to 
the present Hubble rate. We obtain its exact value by 
iteratively solving the equations for inflation and reheat- 
ing. 

We start the calculation for each gravitational wave 
mode from a point where the mode is well inside the 
horizon. The initial condition of the mode is set to oscil- 
late with the amplitude given in Eq. ^. Then we solve 
the evolution of the gravitational waves for each mode 
using Eq. ([2]) and the evolution of the Hubble rate. Af- 
ter the mode has reentered into the horizon and again 
begun to oscillate rapidly, we replace its evolution with 
the WKB solution of Eq. ([7]) since it would be time con- 
suming to treat the oscillations numerically. Note that 
we also use the WKB solution even when the neutrino 
anisotropic stress is taken into account because the effect 
is negligible if fc ^ aH (see Appendix [C| . Finally, we 
compute the spectrum according to Eq. ([5]). 



A. Chaotic inflation 

First, we consider the case where the scalar field has a 
quadratic potential. 



( chaotic inflation: m 56 model ) 



1 




V{(t)) = -m^(J3- 



2 J,2 



(35) 



We used the normalization of the scalar perturbations 
which is derived by combining the constraints from recent 
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FIG. 1: The spectrum of the gravitational wave background 
generated during inflation with a quadratic potential V = 
771(^)^/2, versus the gravitational frequency /o = fc/27r[Hz]. 
The black-line spectrum includes the effects of (7, changes 
and the neutrino anisotropic stress, and the light gray one 
does not. Two lines are shown to compare our numerical re- 
sult with theoretical prediction: the dotted curve represents a 
scale-invariant spectrum and the dashed line is the spectrum 
predicted under the slow-roll approximation, which is plotted 
according to Eq. (|22p with the spectral index nT(k^wot) = 
-1.76 X 10"^ and its running aT(fcpivot) = -3.11 x IQ-". The 
shaded area corresponds to the spectrum due to quantum zero 
point fluctuations (see Sec. IIV A[l . 



observations [42l |. 

^s,prim(fcpivot = 0.05Mpc-i) = (2.45 ± 0.23) x 10"^ 

(36) 
This fixes the value of the scalar field mass to be to = 
1.64 X lO^^GeV. The decay rate is set to be T = 10"^^^ 

Figure [1] shows the result of our numerical calculation. 
In the lower frequency region, we see a well-known char- 
acteristics of the spectrum as described in Sec|TT]- that is 
the change of the frequency dependence of the spectrum 
(r^GW oc fc-^ to fc°) at around 10" ^^Hz, which corre- 
sponds to the horizon scale at matter-radiation equality. 
At frequencies above lO-^^Hz, the spectrum shows some 
interesting structural features, which we will focus on for 
the remainder of this paper. 

Let us first discuss the effect of inflation by paying at- 
tention to the tilt of the spectrum. For comparison with 
previous works, we plotted the Taylor-expanded spec- 
trum derived in the slow-roll regime (Eq. ((22|) ). Note 
that, when comparing with the spectrum of g* = const 
and Hy = 0, our calculation is in good agreement around 
the pivot wavenumber /cpivot = 0.05Mpc- which corre- 
sponds to /o — lO^i^'^Hz. On the other hand, we see 
deviations from the Taylor-expanded spectrum at fre- 
quencies above about 10-^Hz, which includes the fre- 
quency bands of LISA [l^ and DECIGO/BBO [lllll. 
This is because the Taylor expansion around the pivot 
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FIG. 2: Portion of the same spectrum of Fig. [T] focusing 
on the features induced by the change of g, and the neu- 
trino anisotropic stress. Note that Slew is plotted on a linear 
scale. The black spectrum includes both contributions, the 
dark gray one includes only the effect of g, changes, and the 
light gray one does not include either. 



scale (/pivot — 10~^^"'^Hz) is no longer valid. This can be 
avoided if one uses Eq. P^)) instead of Eq. P^ as long as 
the slow-roll condition is satisfied. However, at frequen- 
cies above lO^Hz, where one finds the deviation becomes 
much larger, the slow-roll condition is no longer satisfied. 
This causes a deviation from the correct amplitude, and 
is unavoidable in the slow-roll regime. Therefore, the 
slow-roll approximation overestimates the amplitude of 
the spectrum at very high frequencies. Moreover, the use 
of the Taylor-expanded form of the slow-roll approxima- 
tion, which is used widely in previous works, makes the 
deviation larger and also increases the frequency range 
over which it occurs. Obviously, these deviations are also 
present when the change of 17* is taken into account. 

Other characteristic features in the spectrum are pro- 
duced after the end of inflation. When comparing with 
the spectrum which does not include g* changes and the 
neutrino anisotropic stress, we find these two effects cause 
damping over wide frequencies. While the changes of g^ 
produce damping of higher frequencies (above 10~^^Hz) 
as explained in Sec. IIII Ci the neutrino anisotropic stress 
mainly causes damping at lower frequencies (10~^''Hz to 
10~"'^'^Hz). At frequencies where the two effects overlap, 
the spectrum shows features produced by a combination 
of both the change of g^ and the neutrino anisotropic 
stress. We present a close-up of the overlapped region 
in Fig. [21 Note that we additionally show the spec- 
trum which only includes the effect of g* changes to 
make it clear which effect causes the damping at each 
frequency. It may be seen that the stepwise changes 
of g^ induce a steplike shape in the spectrum and the 
neutrino anisotropic stress only affects the modes which 
reenter the horizon after the time of their decoupling at 
2MeV. The magnitude of the damping due to neutrino 
anisotropic stress may be derived analytically to be 35.6% 



[2IJ . Note that our result does not show the dip and peak 
around /o = 0.8 x lO'i'^Hz (fco = 5 x lO^^'^Hz) as seen 
in the results of Ref. [22] . We find that the feature does 
not arise for the reason given in Appendix D of Ref. 23| , 
but instead is the result of inaccurate treatment of the 
source term. For details, see Appendix [Ul ^ 

The change of the Hubble rate produces features in the 
spectrum around the frequency of the mode that reenters 
the horizon at the time of the change. One can relate the 
frequencies today to the temperature of the Universe at 
the time when the modes reenter the horizon as 1431 



/o = 1.65 X 10"' X — 
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r The 1 
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\9*iThc)] 
100 



Hz. (37) 



For example, substituting the temperature of the quark- 
gluon phase transition (called QGP P.T. in the figures), 
which is assumed to occur at ISOMeV, for The yields 
/o ~ lO^^Hz. Indeed, at around this frequency, we do 
find the large step due to the sudden change of g* [2^ . 
Also, the change of the frequency dependence (I^gw oc fc° 
to k~^) at around /o ~ lO^Hz corresponds to the reheat- 
ing temperature Trh — 2 x lO^^GeV, which is derived 
by substituting F = 10" 



~ lO^^GeV into Eq. ([27|. 
It is notable that the k^'^ dependence above lO^Hz is 
the reheating effect due to the fact that the Hubble rate 
has the same evolution as a matter-dominated Universe 
during reheating. 

One thing we have to mention here is the significant 
rise at the highest frequencies, which is proportional to 
k'^. These are the modes which have not crossed outside 
the horizon during inflation, so that their amplitude de- 
crease as r^GW fx a~^ all the way through keeping the 
gradient of the initial spectrum. Of course they may 
not in truth exist because they does not "freeze in" as 
classical fluctuations. One may regard them not as grav- 
itational waves which have been propagating since the 
epoch of inflation, but as quantum perturbations of the 
space-time metric which are arising "just now" . ^ 



B. LoAV reheating temperature 

Here, we consider different reheating temperatures, re- 
ferring to the upper bound from the gravitino problem, 



We thank Y. Watanabc (private communication) for helping 
to clarify the reason for the difference between our results and 
theirs. 

These modes are considered to be contributing to the cosmologi- 
cal constant and tuned to be very small. However, for reference, 
if they contribute to the total energy density of the Universe, 
the energy density of the gravitational waves exceeds that of the 
scalar field for /o > 10^^ Hz, and the linear analysis becomes 
invalid. 
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FIG. 3: Ttie spectrum for different values of ttie retreating 
temperature for cliaotic inflation. The lightest gray spectrum 
is the same as shown in Fig. [T] i.e. the decay rate is set to be 
r = 10"^m which corresponds to Trh ~ 2 x lO^GeV. The 
dark gray one shows the case of Trh — lO^GeV and the black 
one shows the case of Trh — lO^GeV. 



Trh < 10^ - lO^GeV ^. As we discussed in the previ- 
ous subsection, the reheating temperature determines the 
characteristic frequency of the feature induced by reheat- 
ing, so that different reheating temperatures are expected 
to result in different spectrum shapes. Here, we consider 
the case of Trh = lO^GeV and lO^GeV, which corre- 
sponds to the decay rate F ~ 2 x lO^^^m and 2 x 10~^'^m. 
The calculation is carried out with the same condition as 
that of quadratic chaotic inflation except the value of the 
decay rate. 

Shown in Fig. [3] is a comparison of the spectrum at 
different reheating temperatures. It can be seen that 
each spectrum has the fc~^ dependence due to the re- 
heating effect, although their length are quite different. 
This is because a lower decay rate makes the reheating 
time longer. Note that the edge at the frequency which 
corresponds to the reheating temperature shifts to lower 
frequencies as the reheating temperature decreases. This 
behavior is consistent with the estimation from Eq. (j37p 
that the corresponding frequencies for Trh = lO^GeV 
and lO^GeV are 10~^Hz and 10^ Hz which will be cov- 
ered by DECIGO/BBO. It is an exciting possibility to 
determine the reheating temperature of the Universe by 
these experiments |4^, |4g| . 



C. Application to other inflation models 

We also repeat the calculations for other inflation mod- 
els: A(/)^, new inflation, and hybrid inflation. In each case, 
we adopt the same normalization as in Eq. p6p and set 
the decay rate to be F = 10~^to, where rh is the effective 
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FIG. 4: The spectrum for different inflation potentials. 



mass of the scalar field at the minimum of the potential. 
The results for each model are summarized in Fig. U) 
The values of ut, ax, which are required to draw a line 
of the slow-roll prediction, and r are listed in Table [D 
Also, the comparison between all models is shown in Fig. 

El 



1. A(j!>* model 

Let us consider chaotic inflation with a quartic poten- 
tial. 



1 



v{cj>) = -M\ 



(38) 



where the normalization gives A = 1.78 x lO^^^GeV. A 
notable characteristic of this type of potential is that the 
decay rate is not constant, because the oscillation fre- 
quency of (p depends on the amplitude of (j), i.e. m 
c^/\{(j}), where c is a constant of order 1 [43, Hs 



In the 



case, the Hubble rate behaves as in a 



radiation-dominated Universe when the (j) field oscil- 
lates about the minimum of the potential. Indeed, as 
one can find from the spectrum at the top of Fig. 21 
there is no matter-dominated phase (k~'^ dependence) 
soon after inflation. Therefore, for 0* model, we can- 
not distinguish the stage of reheating and the subse- 
quent radiation-dominated era, so that the uncertainty 
of the decay rate does not affect the resultant spectrum 
at all. Note that the large gap at high frequencies is 
only because of the large deviation from the slow-roll ap- 
proximation, which is also an interesting characteristic 
of this potential. Again, as is seen from the figure, the 
Taylor-expanded spectrum overestimates the spectrum 
at /o > 10-5 Hz. 



S. new inflation 

We now consider the original new inflation model 
which is based on the Coleman- Weinberg potential [43 |. 



y(0)4A^^(ln^-i) + 



16 ' 



(39) 



The oscillation frequency is given hy m — a\X. Here, we 
present the results for two different parameter sets: (1) 
a = IQmpi, A = 2.36 x IQ-i^GeV and (2) a = IQ-'^mpu 
A = 7.27 X IQ-^^GeV. In case (1), we find a very similar 
spectrum to the case of the rn?(fp' potential. By con- 
trast, in case (2), the spectrum is quite different. This 
difference arises from the difference in the point in the 
potential at which the (j) field rolls down. In case (1), the 
scalar field rolled down where the potential is quadratic, 
and this is why the spectrum has almost the same shape 
as that of the rri^cfp' case. 

In case (2), inflation occurs at a lower energy scale, and 
the decay rate is set to be smaller, F = IQ-^m, in order to 



Model 


riT 


(y.T 


r 


m2,^2 


-1.76 X 10^^ 


-3.11 X 10"'' 


0.141 


A</>* 


-3.40 X 10"^ 


-5.78 X 10"'' 


0.272 


New (1) 


-1.16 X 10"^ 


-2.60 X 10-'' 


9.31 X 10-^ 


New (2) 


-6.83 X 10-'** 


-4.77 X 10-'^ 


5.47 X 10-" 


Hybrid 


-5.80 X 10-* 


3.30 X lO-** 


4.64 X 10-^ 



TABLE I: A surninary of the values of the inflationary pa- 
rameters ut, Q.T, and r for each model. These are evaluated 






pivot 



: 0.05Mpc" 



not exceed the Hubble rate during inflation. The remark- 
able features of case (2) are the complete flatness of the 
spectrum and the smallness of its amplitude. The flat- 
ness is due to the small gradient of the potential because 
the initial value of (j) is very near zero. The smallness 
of the amplitude relates to the flatness of the spectrum, 
about which we will explain later (see Sec. IIV C 4p . 



3. hybrid inflation 

In the case of hybrid inflation, the potential is con- 
structed with two scalar fields (j) and a [50, IsH , 



y(0) 



1 

4A 



(Af2 - Ar72)2 + iTO2^2 _^ ig2^V. 



(40) 



We consider the usual case where inflation is driven 
by the (j) field and suddenly ends after the symmetry 
breaking of the a field (waterfall field). To satisfy this 
condition, we adopt the following parameters: A = 1, 
5 = 8 X 10-4, m = 1.5 X 10-^TOp;, M = 1.21 x lO^^GeV. 
In this case, the oscillation at the bottom of the potential 
mainly comes from the (j) field, of which the frequency is 
7710 = gM/\/X, thus we can assume that the radiation 
energy produced during reheating comes from only the 
decay of the </> field. 

The spectrum shown at the bottom of Fig. [H is nearly 
flat and in good agreement with the Taylor-expanded 
spectrum. This is because the quite gentle slope of the 
potential of the (f> fleld results in an almost constant value 
of the Hubble rate as long as the a field triggers the sud- 
den end of inflation. Thus, the spectrum has less fre- 
quency dependence, which leads to less deviation from 
the Taylor expansion. Also, interestingly, the value of 
ax is positive, which means rj > 2e. Both of these are 
the specific features of hybrid inflation. 



4- Comparison 

Finally, in Figl5l we compare all models. In each 
case, the effects of neutrino anisotropic stress and the 
5, changes appear of course in the same way at the same 
position because these shapes are formed well after the 
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const & Ilij = 



2.57 X 10" 



1.26 X 10"^^ 
1.35 X lO"^"^ 



3.39 X 10" 
8.79 X 10" 



1.07 X 10" 

1.21 X lO"^"^ 1.07 X 10" 



2.27 X lO"^*^ 
5.92 X 10"^^ 
8.80 X 10"^'' 



<7* yj- ) iy^ i-l-ij, neutrinos 

g* = const & Ilij = 
Taylor-expanded slow-roll 



3.47 X lO"^'' 
5.01 X lO"^'' 



5.37 X 10"^® 

1.38 X 10"^^ 



3.85 X 10 



-16 
-15 



2.65 X 10" 



1.00 X 10~^^ 6.93 X 10" 



1.58 X 10"^^ 1.27 X 10" 



■15 



1.01 X 10" 



1.24 X lO^^'' 
3.27 X lO"'^'' 
6.99 X 10^^'^ 



4.80 
12.0 
15.2 



Chaotic (m2(/)^ Trh = lO'^GeV) 


Q^\J- ) CZ llij, neutrinos 


2.75 X lO"^'' 


4.90 X 10"^^ 4.01 X lO"^*^ 3.17 x 10"^^ 


9.12 X 10"^** 


4.53 


Chaotic (m^(/>^, Trh = lO^GeV) 


Q'-^\J- ) & iiij, neutrinos 


2.88 X 10"" 


4.90 X 10"^*^ 3.94 X 10"^^ 1.51 x 10"^^ 


7.60 X 10"^'' 


0.20 


Chaotic {\<l)'^) 



3.66 
9.46 
14.2 



New inflation (a = lOrripi) 



5*(T) &nij,nentrinos 1.20x10" 

g. = const & Hij = 1.70 x 10" 

Taylor-expanded slow-roll 



3.89 X 10""'' 
1.00 X 10"" 
1.07 X lO"^"^ 



3.40 X 10"^^ 2.91 X 
8.81 X 10"" 7.55 X 
9.77 X 10"" 8.95 x 



10"" 
10"" 
10"" 



2.09 X 10"" 
5.47 X 10"" 
7.78 X 10"" 



4.36 
11.0 

14.7 



New inflation [a = 10 rn,pi) 



<7* yl } OZ ilij, neutrinos 

g* = const & Ilij = 
Taylor-expanded slow-roll 



7.08 X 10"^^ 
1.00 X 10"^° 



3.55 X 10" 



3.54 X 10" 



3.54 X 10" 



9.02 X 10"^^ 9.01 X 10"^^ 9.01 x 
9.02 X 10"^^ 9.02 X 10"^^ 9.02 x 



10" 
10" 



3.51 X 10"^^ 
8.97 X 10"^^ 
9.02 X 10"^^ 



4.95 X 10"" 
1.27 X 10"" 
2.58 X 10"" 



Hybrid inflation 



g*{T) & nij.ncutrinos 5.75 X 10" 

5* = const & liij = 8.32 x 10" 

Taylor-expanded slow-roll 



1 X 10"" 2.85 X 10"" 2.84 x 

,36 x 



7.41 X 10"" 7.38 X 10"" 7 
7.41 X 10"" 7.38 X 10"" 7.36 



10"" 
10"" 



X 10" 



2.77 X 10"" 
7.18 X 10"" 
7.34 X 10"" 



0.41 
1.02 
1.05 



TABLE II: The density parameter of gravitational wave background r^cw at /o = 10 



,10" 



10"^ 10"\ lO^Hz for several 



inflation models, (1) calculated by fully numerically solving the evolution of gravitational waves including g* change and the 
neutrino isotropic stress, or (2) assuming g, = const and Ilij = 0, or (3) estimated by the Taylor-expanded form of the slow-roll 
approximation using Eq. p2|) with p, = const and Ilij = 0. In addition, the rightmost column shows the signal-to-noise ratios 
for DECIGO, which is aimed at 10"^Hz. 



end of inflation. Needless to say, the differences between 
the models arise due to the difference in how the scalar 
field evolves in each potential during inflation. The be- 
havior of the Hubble rate during reheating is also impor- 
tant as a model dependent factor. Indeed, the reason for 
the unique shape of the Xff)^ model at high frequencies 
is because the Hubble rate evolves differently than the 
other three models. 

Note that the flatter the spectrum, the more the am- 
phtude of large-scale modes is suppressed. This behavior 
can be explained as follows: From Eq. (P5|) . one can 



easily find Vt, 



prim 



= rVs, 



prim ■ 



Since we normalize the 



amplitude of the spectrum using the value of 7^s,prim, 
the smaller scalar tensor ratio r results in a smaller am- 
plitude of the primordial gravitational wave spectrum. 
In addition, r can be rewritten as r ~ —Sut, the tilt of 



the spectrum. This is equivalent to saying that a flatter 
spectrum leads to a smaller amplitude primordial gravi- 
tational wave spectrum. This can be confirmed by com- 
paring the values of r or ut given in Table U with the 
spectrum amplitudes at the pivot scale (see the case of 
g* = const and H^^ = 0), presented in the first column of 
Table HIl The extremely tiny value of r in new inflation 
case (2) is the reason why the amplitude of the spectrum 
in that case is quite small. 

We also provide in Table |ll] the detailed numerical 
values of the spectrum at /o = 10"^, 10"^, 10"^, lO^Hz 
(including the frequency bands of LISA [12| . DE- 
CIGO/BBO [il[il and LIGOH/LCGT [ll,[i5) for sev- 
eral inflation models. For comparison, the values of f^cwj 
calculated assuming g^, = const and H^j = or using the 
Taylor-expanded spectrum predicted under the slow-roll 
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FIG. 5: The spectrum of gravitational wave background is 
shown for several inflation models. The solid curve shows the 
spectrum for the quadratic potential, the dashed curve for the 
quartic potential, the dotted curve for the new inflation po- 
tential with a = lOmpi , and the dot-dashed curve for hybrid 
inflation. 



approximation, Eq. ((22|) . are also given. These clearly 
illustrate the suppression of the precise values compared 
to those derived under these assumptions. For reference, 
we calculate the predicted signal-to-noise ratio (SNR) us- 
ing the noise spectrum of DECIGO, which assumes a 
Fabry-Perot type space interferometer and ten years of 
observing time (Tobs = lOyr) [4l|, [S^]. This is shown in 
the rightmost column in TableHH Note that we find large 
differences in SNR between results from the spectrum we 
have obtained and the Taylor-expanded spectrum. This 
clearly illustrates the reason why it is necessary to eval- 
uate the spectrum amplitude numerically. We find that 
the damping effect due to the changes of g^, decreases the 
SNR dramatically, and also the assumption of the Taylor- 
expanded slow-roll overestimates the SNR in some cases. 
The difference is more prominent for the models with 
steeper potential such as chaotic inflation models. A de- 
tailed analysis of the detectability will be done using the 
results of this paper in our future work. 



V. CONCLUSION 

Gravitational waves generated during inflation are ex- 
pected to be a powerful tool to probe the physics of the 
early Universe. In particular, direct detection of gravi- 
tational waves would provide vital information on their 
high frequency modes. In this paper, we have calculated 
the spectrum of the gravitational wave background over a 
wide range of frequencies (from 10~^^ to lO^Hz) including 
the very high frequencies where we expect to see a devia- 
tion from the slow-roll approximation and the effect of re- 
heating. This has become possible by numerically solving 
the evolution of a scalar field from quantum states during 
inflation until today taking into account the conversion 



of the scalar field energy into radiation. In addition, we 
have also taken into account the contribution of neutrinos 
to the anisotropic stress and the temperature dependence 
of the effective degrees of freedom g* . These two effects 
had already been investigated by Ref. 2^. However, we 
found that a correction is needed for the treatment of the 
source term in their calculation. Moreover their calcula- 
tion is limited in terms of frequency range, so that their 
result is not relevant to the reheating process. Therefore, 
this is the first calculation that correctly includes all the 
factors which affect the shape of the gravitational wave 
background spectrum over a wide range of frequencies. 
This comprehensive treatment of the factors enables us 
to make an accurate estimation of the amplitude of the 
spectrum, which should be invaluable in making predic- 
tions for future experiments. 

Our estimation shows that the Taylor-expanded spec- 
trum of the slow-roll prediction Eq. (j22|) . which is analyt- 
ically obtained and used in many previous works, over- 
estimates the amplitude of the spectrum in some mod- 
els. The deviation between our result and the analyt- 
ical prediction is especially evident at high frequencies 
(/o ^ lO^Hz) because the slow-roll conditions are vi- 
olated towards the end of inflation. It is also notable 
that the deviation is found even in the frequency region 
where the slow-roll conditions are satisfied (around be- 
tween lO^Hz and lO^^Hz in the rn^cjp case). This is be- 
cause the Taylor expansion of the spectrum around the 
CMB scale is no longer a good approximation at these 
frequencies. However, since we do not know the nor- 
malization of the scalar perturbations at much smaller 
scales (< kpc), it is impossible to carry out a Taylor 
expansion about a high frequency point. Therefore the 
deviation is inevitable as long as we analytically evaluate 
the spectrum using the Taylor expansion as in Eq. (j22p . 
We strongly recommend numerical calculation of spectra 
in order to discuss the detectability of the gravitational 
wave background. 

Moreover, we have shown that the spectrum has the 
characteristic features produced after the end of infla- 
tion. The effect of reheating is observed in the spectrum 
as a dip at ~ lO^Hz in the case where the scalar field 
oscillates in a quadraticlike potential. This is because 
the matter-dominated evolution of the Universe during 
reheating causes the frequency dependence of fc~^. This 
feature is not seen in the case of a quartic potential be- 
cause of the difference in the evolution of the Hubble 
rate. In addition, the combination of the effects due to 
changes in g^, and the neutrino anisotropic stress is ob- 
served between 10^^^ — lO^^Hz. It should be noted that 
changes in g* decrease the amplitude by about half an or- 
der of magnitude at higher frequencies, where the above- 
mentioned features of the early Universe notably appear. 
This damping effect becomes larger in case where there 
exist SUSY particles or other exotic particles which cause 
a change of <?* [22] ■ We have confirmed that this effect is 
also important when evaluating the detectability. 

Finally, we would like to conclude with some prospects 
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for future observations. When comparing the results of 
different inflation models, we found that the amplitudes 
of the three models, both cases of chaotic inflation and 
case (1) new inflation, are comparable around the fre- 
quencies covered by future experiments like LISA, DE- 
CIGO/BBO and LIGOII/LCGT. Thus, to distinguish 
these models, it may be necessary to measure the tilt of 
the spectrum. Furthermore, interestingly, the frequencies 
covered by future experiments cover the range in which 
the spectrum shows an edge due to the effect of the low 
temperature reheating. Hopefully, they will provide a 
new insight into understanding of reheating. 

Also, high frequency gravitational wave background 
around lO^Hz, which the new design of detector sug- 
gested in Ref. [Ij] is targeting, will be a probe of re- 
heating process soon after inflation. This is also the fre- 
quency at which we found the connecting point to the 
k^ dependence, beyond which the waves are considered 
to not "freeze in" as classical fluctuations. However, the 
mechanism of " freeze in" is yet to be uncovered. There- 
fore, through future direct detection at around lO^Hz, we 
may obtain information not only on reheating but also 
on classicalization of the quantum fluctuations. 
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APPENDIX A: QUANTIZATION OF TENSOR 
PERTURBATIONS 



To describe the initial condition of gravitational waves 
during inflation, we need a quantum-mechanical treat- 
nrcnt. Here, we do not go into details but give a brief 
description of how to proceed (For details, see Refs. 
[11,12^]). 

The quantization is carried out by introducing the 
rescaled variable ^^(t, k) = ah^/y/WirG which sat- 
isfies the commutation relations with its momentum 



7r>i{T,k)^dvt{T,k)/dT, 






(Al) 



We take the expression for the quantum operator 



^(r, k) = Vk{T, k)at + vUt, fc)a^i, (A2) 



where the creation and annihilation operators, a^ and 

a^arV'5'"'5(3)(k-l), 



a_j^, satisfy the commutation relations 



[«k,«r] 



r-Af -crtl 



0. 



This follow from the Wronskian condition 



i>k- 



dT 



Vk 



<K 

dT 



(A3) 



(A4) 



Then Eq. ([2]) can be recast in terms of conformal-time 
derivatives as 



dT^ 



+ k'- 



Id'^a 



Vk = 0. 



(A5) 



Using the relation between the conformal time and the 
expansion rate during inflation, r ~ —l/{aH) which is 
equivalent with da/dr ~ —aJT, the second derivative 
term can be rewritten as {d'^ a / dr"^) / a ~ 2/r^. When con- 
sidering the short-wavelength limit k/aH — > oo, we can 
neglect the term of 1/r ~ —aH which is much smaller 
than the k term, and obtain the asymptotic solution us- 
ing the normalization of Eq. (|A4[) , 



1 



Vk 



'ikr 



2fc 



(A6) 



This is the solution for the modes well inside the horizon 
during inflation - that is what we set as initial condition 
in our calculation. 



APPENDIX B: THE WKB SOLUTION 

Here we present the derivation of the WKB solution, 
Eq. ([7|). First, we rewrite Eq. ([4]) in terms of conformal 
time r, 



d-'hj 

dr'^ 



2 da 
a dr 



dh^ 

dT 



k^hi = 0. 



One may assume the form of the solution as 



(Bl) 



(B2) 



and consider the case in which the oscillation is very 
rapid compared to the time variation of the amplitude, 
which means dA/dr and d'^A/dr'^ are much smaller than 
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-0.1 




FIG. 6: Time evolution of x'(m = kr) = h^'{u)/h^{0) for the 
mode kr^^dcc = 5.0, which entered the horizon before neutrino 
decouphng. The solid line is the analytical solution, x' = 
—ji{u), and does not include the effect of neutrino anisotropic 
stress. The dotted line, which is almost indistinguishable from 
the analytical solution, represents our numerical result. The 
dashed line is an attempt to reproduce the numerical result 
of Ref. [23 ]. The vertical dashed line represents the time of 
neutrino decoupling. 



dB/dr. By substituting Eq. (lB2| into Eq. JBT]), we ob- 
tain two equations arising respectively from the real part 
and the imaginary part, 



(PA_ 
dr^ 



A 



dB 

rf7 



dAdB 
2 

dr dr 



A 



+ 


'2 da 
a dr 


dA 
dr 


-f 


■s 


B 

T + 


2 da' 
a dr 





k'A 



dB 



0, 

:0. 



(B3) 
(B4) 



Using the assumption that d^A/dr^ ^ dB/dr and con- 
sidering the modes well inside the horizon [da/ dr)/ a <C 
/c, we obtain the solution of B{t) from Eq. (jBSp . 



B{t) ~ ztkr + const. 



(B5) 



The solution of A(t) is derived by substituting this into 
Eq. dH, 



(B6) 



A{t) (X a ^. 

Therefore we obtain the WKB solution, 

C 
/i1(.(t) ~ — exp(±ifcr + a), 
a 

where C and a are arbitrary constants. 



(B7) 



APPENDIX C: ANISOTROPIC STRESS OF 
FREE-STREAMING NEUTRINOS 

Analytically solving the linearized Boltzmann equation 
for neutrinos, we can calculate the spectrum including 
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FIG. 7: Evolution of the source term, i.e. the right-hand side 
of Eq. (IC1[I . The upper figure shows our results and the 
lower shows the reproduced results of Ref. [23 |. The solid, 
dotted and dashed lines are the same as in Fig. [S] and the 
short-dashed line represents the source term. The artificially 
introduced discontinuity in the source term in the lower panel 
has been chosen by trial and error to recover the evolution of 
Ref. [23 ]. Note that the scale of the vertical axis for the 
source term is 10 times larger than that for the others. 



the effect of free-streaming neutrinos by only evolving 
h^. According to Ref. [2l| the evolution equation, Eq. 
^ , can be rewritten in the form of an integro-diffcrential 
equation with using conformal time t. 



d^h^iir) 

dT^ 


'2 da' 
a dr 


dh^ir) 
dr 


X 


'1 da 
a dr 


7: 


dr' 

dec 



j2(k{T~T')) 



-24/, (r) 



B{t-t') 



A2 



dhjir') 
dr 



,(C1) 



where J2 is the spherical Bessel function of order 2. The 
function /j/(t) describes the fractional energy density of 
neutrinos, which is defined using the total energy density 
p and the energy density of neutrinos pv as 



U{r) 



p{r) 



lv,0 



1 -I- a{T)/a, 



(C2) 



eq 



where f^^ ^ fl^/in^ + n^) = 0.40523. 

Neutrinos do not induce an anisotropic stress at ear- 
lier epochs of the Universe, because they are in equi- 
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librium with electrons and photons which are coupled 
to the nonrelativistic baryonic matter whose anisotropic 
stress is negligible. After they decouple and begin free- 
streaming, their anisotropic stress affects the evolution 
of the gravitational waves as viscosity. As seen in Eq. 
(jCip . the source term depends on the Hubble rate, which 
means the anisotropic stress does not affect the mode 
when it is deep inside the horizon (k 3> aH). Thus, 
the damping effect does not arise at the higher-frequency 
modes which were already inside the horizon at the time 
of neutrino decoupling. Also, once the Universe be- 
comes matter-dominated, the contribution of neutrinos 
to the total energy density of the Universe becomes suf- 
ficiently small (/jy(T) <C 1), thus the damping effect does 
not appear in the modes which enter the horizon dur- 
ing the matter-dominated era. Therefore, the damping 
effect of the neutrino anisotropic stress is seen between 
the modes which correspond to the time of neutrino de- 
coupling (/o ^ 10^ ^"^Hz) and radiation-matter equality 
(/o ~ 10-i^Hz). ^ 

Here, we explain the reason why our spectrum does not 
have the dip and peak which are seen in the spectrum ob- 
tained by Watanabe and Komatsu ^^. In their paper, a 
small ripple is seen around the frequency corresponding 
to the time of neutrino decoupling, /o = 0.8 x 10^^"Hz 
{ko = 5 X 10~"'^°Hz), and is interpreted as a result of as- 
suming instantaneous decoupling. As seen in Fig. [21 our 
result also shows a ripple at around the same frequency, 
however it does not show a large dip and peak as ob- 



served in their spectrum (see Fig. 5 of Ref. [23 | ). The 
detailed reason for the peak is explained in Ref. |22l | us- 
ing a plot of the time evolution of the amplified mode 
(see Fig. 10 of Ref. [22]). Their suggestion is that the 
mode gains energy from neutrinos when the derivative 
of the mode, x'(u = kr) = h^{u)/h^{0), is positive, so 
that it is amplified in the case where the mode enters the 
horizon when x' > 0. However, we have reproduced the 
same behavior by changing the source term, and found 
that the amplification does not arise for the reason sug- 
gested in their paper. 

The reproduced behavior of the derivative of the mode 
is plotted in Fig. [61 which presents the same quantities 
as Fig. 10 of Ref. [22|. We also plot the result of our cal- 
culation in the same figure. Note that the amplification 
which is seen in the reproduced mode does not appear 
in our result. This difference is explained by the behav- 
ior of the source term, which is presented in Fig. [T] As 
seen in Fig. [71 the source term has an unnatural discon- 
tinuous jump in the case where the mode is amplified. 
This behavior can be reproduced when we calculate the 
source term by replacing the lower limit of the integral 
in Eq. (jCip with zero and switch the source term to af- 
fect the evolution of gravitational waves at the time of 
neutrino decoupling. Therefore the dip and peak seen in 
Ref. |22l | are due to inaccuracy in computing the source 
term. The validity of our calculation has been confirmed 
by the author of Ref. 2^, Y. Watanabe, via private 
communication. 



[1] A. H. Guth, Phys. Rev. D 23, 347 (1981). 

[2] A. Albrecht and P. J. Steinhardt, Phys. Rev. Lett. 48, [15] 

1220 (1982). 
[3] A. Linde, Phys. Lett. B 108, 389 (1982). 
[4] A. H. Guth and S. Y. Pi, Phys. Rev. Lett. 49, 1110 

(1982). [16] 

[5] A. D. Linde, Phys. Lett. B 116, 335 (1982). [17] 

[6] J. M. Bardeen, P. J. Steinhardt, and M. S. Turner, Phys. 

Rev. D 28, 679 (1983). 
[7] B. Allen, Phys. Rev. D 37, 2078 (1988). 

[8] V. Sahni, Phys. Rev. D 42, 453 (1990). [18] 

[9] L. P. Grishchuk and Y. V. Sidorov, Phys. Rev. D 42, 3413 

(1990); L. P. Grishchuk and M. Solokhin, Phys. Rev. D [19] 
43, 2566 (1991). 
[10] U. Seljak and M. Zaldarriaga, Phys. Rev. Lett. 78, 2054 [20] 
(1997); M. Zaldarriaga and U. Seljak, Phys. Rev. D 55, [21] 

1830 (1997). [22] 

[11] M. Kamionkowski, A. Kosowsky, and A. Stebbins, Phys. 

Rev. Lett. 78, 2058 (1997); Phys. Rev. D 55, 7368 (1997). [23] 

[12] P. Bender et a/., "LISA Pre-Phase A Report" (second 

edition) (1998). [24] 

[13] N. Seto, S. Kawamura, and T. Nakamura, Phys. Rev. [25] 
Lett. 87, 221103 (2001); S. Kawamura et al, Class. 
Quant. Grav. 23, S125 (2006). [26] 

[14] S. Phinney et al, NASA Mission Concept Study, 

http://universe.nasa.gov/prograni/bbo.htnil ; G. M. [27] 
Harry, P. Fritschel, D. A. Shaddock, W. Folkner and E. S. 
Phinney, Class. Quant. Grav. 23, 4887 (2006) [Erratum- [28] 



ibid. 23, 7361 (2006)]. 

E. Gustafson, D. Shoemaker, K. A. Strain, and 

R. Weiss, LIGO Document No. T990080-00-D 

(www.ligo.caltech.edu/docs/T/T990080-00.pdf); P. 

Fritschel, Proc. SPIE, 4856, 282 (2003). 

K. Kuroda et ai, Class. Quantum Grav. 23 S215 (2006). 

A. Nishizawa et al, Phys. Rev. D 77, 022002 (2008); A. 
Nishizawa et al., Classical Quantum Gravity 25, 225011 
(2008); T. Akutsu et al., Phys. Rev. Lett. 101, 101101 
(2008). 

D. Baumann, P. Steinhardt, K. Takahashi and K. Ichiki, 
Phys. Rev. D 76, 084019 (2007). 

C. Carbone, C. Baccigalupi and S. Matarrese, Phys. Rev. 
D 73, 063503 (2006). 

D. J. Schwarz, Mod. Phys. Lett. A 13, 2771 (1998). 
S. Weinberg, Phys. Rev. D 69, 023503 (2004). 

Y. Watanabe and E. Komatsu, Phys. Rev. D 73, 123515 

(2006). 

L. A. Boyle and P. J. Steinhardt, Phys. Rev. D 77, 063504 

(2008). 

M. S. Turner, Phys. Rev. D 55, R435 (1997). 

T. L. Smith, M. Kamionkowski and A. Cooray, Phys. 

Rev. D 73, 023504 (2006); 78, 083525 (2008). 

B. C. Friedman, A. Cooray and A. Melchiorri, Phys. Rev. 
D 74, 123509 (2006). 

J. R. Pritchard and M. Kamionkowski, Annals Phys. 

318, 2 (2005). 

M. Maggiore, Phys. Rep. 331 (2000) 283-367. 



15 



[29] J. E. Lidsey et al, Rev. Mod. Phys. 69, 373 (1997). [41] 

[30] A. Kosowsky and M. S. Turner, Phys. Rev. D 52, R1739 

(1995). [42 

[31] A. D. Dolgov and A. D. Linde, Phys. Lett. B 116, 329 [43 

(1982). 
[32] L. F. Abbott, E. Farhi and M. B. Wise, Phys. Lett. B [44 

117, 29 (1982). 
[33] E. W. Kolb and M. S. Turner, The Early Universe (West- [45 

view Press, Boulder, CO, 1990). 
[34] L. Kofman, A. Linde and A. A. Starobinsky, Phys. Rev. [46 

D 56, 3258 (1997). [47 

[35] J. H. Traschen and R. H. Brandenberger, Phys. Rev. D 

42, 2491 (1990). [48 

[36] L. Kofman, A. Linde and A. A. Starobinsky, Phys. Rev. 

Lett. 76, 1011 (1996). [49 

[37] R. Easther, J. T. Giblin and E. A. Lim, Phys. Rev. Lett. 

99, 221301 (2007). [50 

[38] J. Garcia-Bellido and D. G. Figueroa, Phys. Rev. Lett. [51 

98, 061302 (2007); J. Garcia-BeUido, D. G. Figueroa and 

A. Sastre, Phys. Rev. D 77, 043517 (2008). [52] 

[39] J. F. Dufaux, A. Bergman, G. Felder, L. Kofman, and J. 

P. Uzan, Phys. Rev. D 76, 123517 (2007). 
[40] D. N. Spergel et al, Astrophys. J. Suppl. 170, 377 (2007). 



T. Chiba, Y. Himemoto, M. Yamaguchi and J. 

Yokoyama, Phys. Rev. D 76, 043516 (2007). 

U. Seljak et al, Phys. Rev. D 71, 103515 (2005). 

M. Kamionkowski, A. Kosowsky, and M.S. Turner, Phys. 

Rev. D 49, 2837 (1994). 

M. Kawasaki, K. Kohri and T. Moroi, Phys. Lett. B 625, 

7 (2005). 

K. Nakayama, S. Saito, Y. Suwa, and J. Yokoyama, J. 

Cosmol. Astropart. Phys. 06 (2008) 020. 

S. Kuroyanagi et al, in preparation. 

M. Morikawa and M. Sasaki, Prog. Theor. Phys. 72, 782 

(1984). 

Y. Shtanov, J. Traschen and R. Brandenberger, Phys. 

Rev. D 51, 5438 (1995). 

S. Coleman, and E. Weinberg, Phys. Rev. D 7, 1888 

(1973). 

A. Linde, Phys. Rev. D 49, 748 (1994). 

J. Garcia-Bellido, and A. Linde, Phys. Rev. D 57, 6075 

(1998). 

H. Kudoh, A. Taruya, T. Hiramatsu and Y. Himemoto, 

Phys. Rev. D 73, 064006 (2006). 



